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We discuss how to exploit stochastic resonance with the methods of statistical theory of decisions. To do 
so, we evaluate two detection strategies: escape time analysis and strobing. For a standard quartic bistable 
system with a periodic drive and disturbed by noise, we show that the detection strategies and the physics of the 
double well are connected, inasmuch as one (the strobing strategy) is based on synchronization, while the other 
(escape time analysis) is determined by the possibility to accumulate energy in the oscillations. The analysis 
of the escape times best performs at the frequency of the geometric resonance, while strobing shows a peak 
of the performances at a special noise level predicted by the stochastic resonance theory. We surmise that the 
detection properties of the quartic potential are generic for overdamped and underdamped systems, in that the 
physical nature of resonance decides the competition (in terms of performances) between different detection 
strategies. 
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I. INTRODUCTION 

It is well known that Stochastic Resonance (SR) can be 
exploited, under suitable circumstances, to improve detec¬ 
tion However, the knowledge that the system response 

could possibly be enhanced at special noise level is not suf¬ 
ficient to identify the best method to be employed for signal 
detection. In fact, several indexes have been proposed to quan¬ 
tify the response of a system when SR occurs: e.g., Fourier 
spectrum analysis m, m, synchronization measures through 
strobed data amplification of a reaction coordinate 0, 
escape times processing (Ml. Therefore, it would be inter¬ 
esting to know in advance which method will actually work 
better. Put it in another way, it is interesting to know if a link 
exists between the physics characteristics of the system and 
the performance behavior of feasible detection strategies. 

Let us recall some basic concepts of SR from the point of 
view of statistical decision theory. To be specific, we consider 
a signal S{t) that is a mixture of a periodic drive of amplitude 
a, (angular) frequency w, initial (unknown) phase (pQ, and a 
random (uncorrelated Gaussian) perturbations f: 

S{t) = a sin {ujt 4-(po) + ^{t) (1) 

The conventional wisdom is that Stochastic Resonance (SR) 
occurs when the response of a nonlinear system to the signal 
S{t) can be enhanced at a special noise level. In the standard 
analysis of SR |l2l|3][ll, the response is evidenced by the be¬ 
havior of a component of the outgoing Fourier spectrum, that 
is maximized at a certain noise intensity. The upsurge of the 
Fourier component makes SR appealing for signal detection, 
inasmuch as it is conceivable to exploit the Fourier analysis to 
reveal the presence of an injected deterministic signal M- In 
this scheme, one hopes that a particular combination of noise 
and periodic drive makes it easier to detect the signal, for in¬ 
stance because a threshold is only passed with the help of the 


noise m (although it has been proposed to exploit stochastic 
resonance also for suprathreshold deterministic signals ifTTIl ). 
However, the very concept that noise could be beneficial is 
counterintuitive, if not controversial nnum, in that it sounds 
against good sense that more disturbances can, in the end, im¬ 
prove the detection of a signal. This Gordian knot has been 
cut by noticing that SR only helps detection in suboptimal 
systems ifT^ [fSll . and suboptimal threshold choices ifThll . In 
less general terms, let us suppose we want to decide about the 
presence of a periodic drive; then decision theory proves that 
the Neyman-Pearson scheme is optimized by the Likelihood 
Ratio Test (LRT). To decide about the presence of a periodic 
forcing corrupted by white Gaussian noise LRT amounts to 
the scalar product of the signal S{t) and the mask, the peri¬ 
odic component sin {ujf). This is the optimal detection strat¬ 
egy - the matched filter applied to the input signal S{t) iflTl - 
that cannot be improved adding noise. In practice, one is often 
frustrated in the application of LRT optimal technique and em¬ 
ploys suboptimal strategies where SR based signal detection 
enhancement can genuinely occur. As preeminent examples 
where the optimal strategy is not feasible, we can mention the 
cases where the full trajectory S{t) is difficult to retrieve (as in 
very sensitive Fabry-Perot pendulums for gravitational waves 
detection El [HI), or it is just not available for measurements 
(as in Josephson junctions, for the quantum mechanical nature 
of the dynamical variable ll20H22ll ). In some other cases 
the recorded signal is far too long to be analyzed with optimal 
methods, and the LRT is de facto not applicable. In still other 
cases a nonlinearity transformation in the system allows for 
the occurrence of typical SR pattern El- To visualize the dif¬ 
ficulty, we can imagine to tackle the original problem where 
SR arose; the study of climate changes, based on geological 
evidence of the alternate of dry and cold periods ll24l . Can we 
access the instantaneous temperature of the Earth, that is the 
signal S{f)l Unfortunately, we can just estimate the passage 
from an ice age to a dry age, i.e. the escape time from stable 
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FIG. 1: Sketch of the escape process. The separatrix for the poten¬ 
tial Eq. 0 is at a: = 0. The limits 9 and —9 are discussed in the 
Appendix. 


climate configurations. Even if better techniques were avail¬ 
able to estimate the yearly Earth temperature, say from the 
maximum extension of continental ice sheets, this amounts 
to measure each year the temperature in the coldest day. Put 
another way, one could only observe the so-called strobed dy¬ 
namics obtained illuminating the system at some time inter¬ 
vals. The two sampled dynamics we have just mentioned - 
Escape Times (ET) from dry to cold periods and Strobed Dy¬ 
namics (SD) at prescribed time intervals - are suboptimal, in 
that the optimal LRT strategy, equivalent to the matched filter, 
requires to exploit the whole trajectory (i.e., the full informa¬ 
tion content, or the instantaneous temperature) of the input 
signal. Only for the suboptimal strategies, based on the re¬ 
duced data (e.g. the escape times or the strobed dynamics), 
SR can occur. Indeed it has been shown that; i) noise can be 
used to enhance signal detection through the analysis of ETs 
in the first order standard bistable potential fSl and in a sec¬ 
ond order washboard potential ll25l ii) with the appropriated 
choice of noise intensity SD exhibits good detection perfor¬ 
mances, « ‘idB below the optimum, when strobing occurs at 
the forcing period 27r/a; @|. 

Thus the alternatives to the matched filter are practically 
attractive, and could possibly exhibit a bona fide enhancement 
when noise increases. If one backs down the optimal matched 
filter, physical intuition suggests to seek for best performances 
in the parameter region where SR occurs. Put it another way, 
if an opportunity to improve the analysis by adding noise is to 
exist, one guesses that the resonant condition of SR is the first 
place where to look for such a chance. 

Still, an open problem can be posed; which technique, 
among the many suboptimal ones, best performs in a spe¬ 
cific physical system? The objective of the present work is to 
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FIG. 2; Performances of the detection for the overdamped system as 
measured by the KC index dxc for the ET and SD strategies. The 
figure demonstrates a decrease of the performances of the escape, 
without resonance or peak, and independent of the frequency of the 
periodic drive cu, in the considered noise interval (0.1 < D < 0.4). 
The peak of the strobed technique decreases when the frequency in¬ 
creases, see also Fig. 


systematically characterize the two detection performances of 
two reduced data - the escape times and the strobed dynamics 
- for the overdamped and underdamped prototypal system of 
a quartic potential. The analysis of the performances demon¬ 
strates that the two techniques display different properties, and 
that the best performing strategy depends upon the underlying 
dynamical nature of the resonance. We remark that we have 
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selected two particular strategies which are suitable when the 
full trajectory is not easily available. Other suboptimalities 
can be devised, for example threshold detectors, similar to a 
neuron, that are most suited in biological applications iniiia. 

The paper is organized as follows: in Sect. |^we outline 
the model equations. In Sect. |I^we establish the suboptimal 
detection strategies for ET and SD, that are applied and eval¬ 
uated for the prototypal quartic potential in Sect. |IV] Sect. |V] 
concludes with the physical consequences of the above analy¬ 
sis. 


II. MODELS FOR BISTABLE SYSTEMS 
A. First Order Model 

Let us consider the signal Q applied to a prototypal quartic 
bistable potential; in normalized units the system is governed 
by the following stochastic differential equation a: 

— X = asinltuf -f ipo) + ^(t) = S{t). (2) 

at 

Eq. a is called overdamped because is the high friction limit 
of a nonlinear oscillator driven by a mixture of deterministic 
signal and noise, as will be discussed in the next Subsection. 
In Eq. © a random terms appears, ^(f), to model an addi¬ 
tive noise with autocorrelation function of intensity D, viz 
= 2D5{t — t'), corrupting the external sinusoidal 
drive. Thermal noise has been supposed uncorrelated with ex¬ 
ternal noise and included in the overall noise D. Eor practical 
applications the physical nature of the noise sources (either in¬ 
trinsic thermal noise or external noise) does affects the mathe¬ 
matical treatment of the detection strategies (if the sources are 
uncorrelated). We focus on the simple case of additive noise 
(that is a paradigm in signal processing); however we expect 
that the results can be indicative of the behavior also for other 
noise sources and systems, as non-Gaussian Ezlllll, colored 
noise ll29l[^ in inertial and overdamped systems OTIl . com¬ 
binations of time delay and non-Gaussian noise Il32l - l34l and 
with different mechanisms llT5l[3^ . 

The SR frequency is (approximately) the following: 

(t)o ~ —. (3) 

Here (r) is the average of the ETs when the noise intensity 
is D (angular brackets denote time average, the subscript 0 
indicates the absence of the external drive). 

The quartic bistable potential of Eig. [T] is associated with 
Eq. 0: 


U[x) = 



( 4 ) 


If the periodic deterministic signal is absent (a = 0), for 
low noise (D « AU) the escapes occur at a rate llJTll 


(^) ^ exp 



( 5 ) 


(a) 



(b) 
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FIG. 3: Comparison of different indicators of SR, obtained by sim¬ 
ulations of Eq. 0. (a) The average amplitude of the displacement, 
Eq.^, that decreases with the external frequency, (b) Distortion 
of the ETs distribution, Eq.j28[>, compared to the displacement for 
the middle frequency, a; = 0.2 (the data have ben smoothed to evi¬ 
denced the peak maximum and the minimum). This is the standard 
SR, as the noise at which the peak occurs increases with the external 
frequency. Parameters are the same as in Eig. 
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(a) 



(b) 



FIG. 4: Comparison of different indicators of SR, obtained by simulations of Eq.lj^: Distortion of the ETs distribution, Eq.l|28^, and average 
amplitude of the displacement, Eq.j27|>. (a) Low frequency behavior, w ^ 1. The indexes exhibit a maximum in correspondence of the best 
performances, see next Eig. (b) Geometric resonance frequency behavior, cu ~ 1. The indexes exhibit a monotonic behavior typical of the 
regions without stochastic resonance. 


where tx is the Kramers prefactor 1221, that also gives the 
order of magnitude of the relaxation time of local equilibrium 
ll38l . The normalized barrier height is AU = 1/4 for the 
normalized potential ( 0 . Thus, when the external drive is 
absent, Eq .0 represents the expected average. 

The average escape time (r) can be directly measured in 
experiments Ul |5l, and in many cases it is the only physical 
quantity promptly available lfT8ll20ll2Tl . In fact in some phys¬ 
ical systems as Fabry-Perot pendulums or Josephson junctions 
it is difficult (or impossible) to follow the entire tra¬ 
jectory of the system, while it is possible to detect the passage 
across a separatrix (and therefore the sequence is more rig¬ 
orously defined as a passage time llJTll L In the presence of 
the drive, the average passage time is altered by the sinusoidal 
term. Loosely speaking, one could therefore distinguish if the 
deterministic signal is present {a ^ 0 ) or absent {a = 0 ) 
by means of the escape rate, as the two cases correspond to 
(t) / (t)o or (r) = (r)o, respectively (a rigorous formula¬ 
tion of the problem will be discussed in Sect. III). SR leads to 
a more effective signal detection if an increase of the noise D 
facilitates the distinction between the two responses (with and 
without external drive). 


B. Second Order Model 

The bistable model Eq.0 can be straightforwardly gener¬ 
alized introducing inertia and damping ( 7 ), that lead to the 


following model; 


^ + 7 ^ -x + x^ = asm{ujt + ipo) + ^{t) = S{t). ( 6 ) 


Equation 0 is a second order stochastic differential equa¬ 
tion of a nonlinear bistable oscillator in a potential given by 
Eq.0. For 7 < 1 this amounts to an underdamped or iner¬ 
tial system. The normalizations of Eq. 0 are not the same as 
per Eq.0, thus the two systems cannot be directly compared; 
however, we will focus on some features that are not affected 
by the different normalizations. 

The analytical treatment of the stochastic differential Eq.0 
is much more difficult than the overdamped case, Eq.0, and 
a general solution valid for all dissipation values has not been 
found, yet 1^ . Even in the absence of the periodic drive, 
the unperturbed escape rate is only obtained with some ap¬ 
proximations, therefore the theoretical understanding of sec¬ 
ond order SR is less complete than in the case of first order 
system. At a simple level, however, the interpretation of SR 
for signal detection is the same: subject to a sinusoidal excita¬ 
tion the escape rate moves away from the unperturbed value, 
and such displacement is enhanced for special values of the 
random intensity 1291 . 
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III. DETECTION STRATEGIES 


To fully define the detection strategies, it is usual to formal¬ 
ize the problem as a binary hypothesis test: 

0 : sinusoidal drive is absent 
Til : sinusoidal drive is present 

For this decision problem two different error probabilities 
arise: 

• Iht false alarm probability Pf, also called Type I error 
probability, i.e. the probability to decide for the hypoth¬ 
esis Hi when Hq is true; 

• the miss probability Pm, also called Type II error prob¬ 
ability, i.e. the probability to decide for the hypothesis 
Ho when Hi is true. 

We hrst consider the case of a bistable system where the 
normalized potential U(x) and the normalized noise standard 
deviation D are known and do not depend on the particular 
hypothesis Ho or Hi in force. We also assume that the signal 
parameters (i.e. a, oj, and (po) are known under Hi hypoth¬ 
esis. In this setup the Neyman-Pearson lemma Oil identifies 
the LRT as the optimal detection strategy, for it minimizes, 
among all possible tests, the miss probability Pm at a fixed 
false alarm level Pf. Thus, for a fixed mean observation time 
(T), if we collect M observations y = supposed 

to be independent and identically distributed, the test statistics 
can be written as: 


friUklHi) 

to Uo 


(7) 


where /y(-|'Ho,i) are the Probability Density Functions 
(henceforth PDF) of the observation y under the hypothe¬ 
sis Ho,i, while is a suitable threshold that returns a fixed 
false alarm level. To simplify the computation of the statis¬ 
tics it is useful to compare the normalized natural loga¬ 
rithm of the likelihood ratio (henceforth LLR) with a thresh¬ 
old C = log(C'): 


M-1 

A(y) = XI 


/y(t/fc|7fi) . 

.fY{yk\Ho)_ S 

7X0 


( 8 ) 


The advantage of Eq. 0 is that the statistic A(r) can be 
computed as the sum of the random samples C — {C,k}^SQ , 
that are obtained from the observation via the optimal (in the 
Neyman-Pearson sense), transformation 


Ck 


log 


fY{yk\Hi) 

jY{yk\Ho) 


(9) 


that contains the information of both PDFs, /y(-|’Hi) and 
/y('|7fo)- 


To further simplify the analysis we employ the Kumar- 
Carroll (KC) index cIkc HO): 

tiirc(A) = I (A)i - (A)o ^^^ 

(a2(A)i+a2(A)o) 

where (A)i,(A)o are the estimated average LLR over a pre¬ 
scribed (T), with and without drive, respectively. We have 
also denoted with (T(A)i,tT(A)o the corresponding estimated 
standard deviations. Index dxc is one possibility among 
many MM to differentiate between Hi and Ho and to sum¬ 
marize the detection performances in a single value, and there¬ 
fore often used to quantify SR ifThll . 

The use of the dKC index relies on the following consider¬ 
ations. The performances of a detector are often represented 
by using the Operating Characteristics (OC) of the detector, 
e.g. the curve that describes the behavior of Pm as a func¬ 
tion of Pf. A first simplification is to compute a particu¬ 
lar point of the OC, in which the two error probabilities are 
equal {Pf = Pm), and to call this value error probability P^. 
Moreover, under suitable hypothesis and for large sample size 
M, LLR is asymptotically normally distributed due to central 
limit theorem (see iS). Thus, after simple algebra, we can 
write: 


Pe 


-erfc 

2 



A(A)2d^c(A)\ 

4 2 V 2 J’ 


( 11 ) 


where A(A) = 2 | ct(A)o - cr(A)i | / | a{A)o + cr(A)i |. 
Equation fn) is a decreasing function of A, and therefore ne¬ 
glecting the difference among standard deviations (if it exists) 
it is possible to retrieve an upper bound of P^ that is only 
function of dxc, i-e- 



erfc 


(dKc{A)\ 

\ 2 V 2 ) 


( 12 ) 


The inequality underlines the heuristic character of the 
KC index as an indicator of the detector performance. We also 
note that dKc{A) can be derived from the KC index computed 


on the single Ck, k S [0, M— 1], as dKc{A) = \/MdKc{P)- 
To summarize: the presence of a coherent drive can be de¬ 
tected collecting a suitable numbers of observations, and the 
performances of the detection can be evaluated using the index 
given by Eq.dOi. In the following subsections we compute 
the LLR statistic defined in Eq.0 for two signihcant types of 
detection strategies based on the escape times and the strobed 
dynamics. 


A. Escape time sequence-based strategy 

Let us now analyze detection strategies based on escape 
time sequences M- More specihcally, we examine the run¬ 
ning state of the bistable oscillator, i.e. we retrieve the escape 
sequence (driven or not by an external deterministic signal) 
letting the system to freely evolve. A single escape time 
is, loosely speaking, the time to pass from a basin to the other 
(see the Appendix for the details). 
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In the described process many time scales occur ll45l . The 
particle initial position evolves and relaxes in the proximity of 
the stable points a; = ±1 on a time scale that is governed by 
the inverse of the friction ~ I/ 7 . A second time scale is the 
principal escape rate, as it describes the transition from one 
side of the barrier to the other. Finally, a third time scale, the 
response relaxation time, corresponds to the time necessary 
to achieve statistical equilibrium, that we refer to as running 
state. The running state corresponds to the system free evolu¬ 
tion, when the influence of the initial conditions is negligible. 
In the following of the paper we consider the sequence of es¬ 
cape times from a basin to the other in statistical equilibrium. 
Details of the numerical procedure to retrieve such sequence 
taking into account the other time scales is discussed in the 
Appendix. The procedure generates an ET sequence {Ti}'^i 
that is the starting point of our analysis. 

Other choices are possible ll^ . for instance the system can 
be prepared in an initial (e.g. x = —1) state with a known 
signal phase tpo, and the escape is defined as the shortest time 
to reach the separatrix. Once the separatrix is passed, to pre¬ 
pare the prescribed initial state requires an additional time to 
set the initial condition at each passage with the required sig¬ 
nal phase. The additional time should be included in a careful 
analysis of the performances for signal detection. Moreover, 
the restart process introduces a further complication in the ex¬ 
perimental setup. Thus, for sake of simplicity, we only con¬ 
sider the free running dynamics. 

With the free running procedure one can approximate the 
probability density of the escape times, The distri¬ 

bution depends on the parameters of the deterministic signal 
and of the noise ll46l . and it can be computed for each of the 
two hypothesis 'Ho and Hi to estimate /r('|Ho) and /t-(-|Hi) 
using a non-parametric statistical technique such as the Kernel 
Density Estimation (KDE) ia, that employs a large number 
of training samples (> 10 ® trials) to have stable results (see 
also ED). A finite sequence t = {Tk\^^Q of M retrieved 
ETs (that plays the role of the observations y) is employed 
to decide if a coherent drive is embedded in the perturbation 
[equationwise, to decide if a 7 ^ 0 in Eq.Q]. The decision 
statistic results in the following test: 


M-l 

Hf) = 




/r('rfclHi) 

/r('rfclHo) 


Hi 

< c. 

Ho 


(13) 


Thus, it is possible compute the KC index for the ETs: 

dfcc = dKcW- (14) 

The choice of M is related to the mean escape time without 
drive ( t)o and to corresponding mean observation time (T)o 
by the approximated relation 

{T)o = M{t)o. (15) 


construct a sequence x{tk), to decide if a drive is present. 
The analysis of a strategy based on sampling at constant time 
intervals leads to the following asymptotic result: for very 
short strobing time. At << 2x1^}, the detection amounts to 
the scalar product (or the Eourier analysis, equivalent to the 
matched filter), i.e. it is equivalent to the optimal solution a. 
It is therefore interesting to investigate strobing at finite time 
intervals. It is also known from the theory of SR El that for 
At = tt/oj (sampling at half the deterministic signal period) 
there exists a suitable phase ip where the probability distribu¬ 
tion of the sequence x{tk) peaks at cc = 1 for even k, and at 
X = —1 for odd k. Eormally, computing N = (T)/Tg, where 
(as usual) (T) is the mean observation time and Tg = 27r/w, 
it is possible to define 2N sampling times: 

\ / 2k + 1 \ 

tk = -i — - — n-ipj, fc = 0,1,2, ...,2At - 1 (16) 

to obtain the sequence: 

Xk = i-l)^x(tk). (17) 

In a perfectly synchronized state the system crosses the sep¬ 
aratrix each half period; if this is the case, the system exactly 
jumps from —1 to 1, and back from to 1 to —1. Thus the se¬ 
quence 0 only contains positive elements. On the contrary, 
for a purely random sequence that casually moves in the phase 
space, the sequence ( [T7] i contains as many positive as negative 
symbols. Put in another way, using the Heaviside 0 function 
it is possible to obtain the observations 


Sk = 9(xk). (18) 

These observations are a Bernoulli random variable with pa¬ 
rameter = Prob(sfc = 1). Thus, the LRT decision statis¬ 
tics is given by counting the number of plus signs, i.e. 

2N-1 

N+= Y. ( 19 ) 

k=0 

where the number of observation is M = 2N. As usual, the 
statistics N+ can be compared to a threshold E to decide if a 
signal is present, in that one expects N+ = N for a purely 
random signal, and = 2N for a sinusoidal drive with¬ 
out noise. The above statistics is referred to as strobed sign¬ 
counting stochastic resonance detector (SSC-SRD) ||4|. 

The SSC-SRD performance is described by the false-alarm 
and false dismissal probabilities El: 

Pf = Prob {N+ > T\Ho} = /i/ 2 (r + 1,27V - E), (20) 


B. Strobed Detection Strategy 

Strobing amounts to analyze the solution of either Eq. •El 
or Eq. at constant time intervals t = tk = kAt, i.e. to 


= Prob {N+ < r|Hi} = 1 - (E +1, 2iV - E), (21) 

where Ip{x, y) is the regularized incomplete Beta function of 
order p. The related OCs are typically worse by « 3dB than 
the corresponding OCs of the matched filter El- 
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FIG. 5: Performances of the detection for the second order system, as measured by the KC index (Ikc, for the ET and SD strategies at the SR 
frequency (a) and geometric resonance (b). SR occurs at D ~ 0.025, as per Eq. The behavior in (a) is the analogous of SR for first order 
systems, Fig.[^ The best performances are obtained from the strobing technique in (a), and from the ET in (b). 
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FIG. 6: Performances of the ET (a) and SD (b) strategies of the un¬ 
derdamped system at low external frequencies uj given by Eq.j^. 
The position of SR frequency as a function of the noise intensity is 
reported in the inset. The performances of the SD in (b) are weakly 
dependent upon the drive frequency, and therefore it is not possible 
to clearly identify the behavior of the peak. 


EIG. 7: Detailed examination of (Ikc for the ET (a) and SD (b) 
strategies at the geometric resonance for the underdamped system 
as a function of the deterministic signal frequency. The dashed 
line in (a) denotes the asymptotic performances of the strobed tech¬ 
nique, as displayed in (b). In (a) the escape strategy performances 
improve with a suitable choice of the frequency of the deterministic 
signal. The best performances are achieved at the resonance of the 
system.The performances of the SD detection strategy saturate when 
noise increases (b). 
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FIG. 8: Underdamped escape times strategy for various external fre¬ 
quencies uj. The best performances are observed at a; = 0.8 in (a) 
1501 . in fact the peak performances increase in (a) (tu < 0.8), and 
decrease in (b) (u) > 0.8). This is in contrast with the monotonic 
behavior of the relation the temperature at which the peak oc¬ 
curs for different frequencies first increases (a) and then decreases 
(b). This clearly defines a standard geometric resonance with a fa¬ 
vorite escape time. The first order analogue, Fig|^ is monotone in 
the external frequency to. 


Also in this case, for high values of 2N it is possible to 
approximate the binomial distribution via a normal one (by 
using the well-known de Moivre-Laplace theorem 1481 ) and 
to introduce a suitable KC index for the strobed statistics, i.e. 

dfcc = dKciN+). (22) 

It is noticeable that in this case there is a simple expression 
for dffc', i.e. 




2F+ - 1 

+ 2P+(1 - P+) 


(23) 



FIG. 9: (a) The frequency dependence of the normalized escape time 
of the underdamped system for two values of the dissipation param¬ 
eter ({r)o is the average escape time in the presence of a purely 
random signal). The dip occurs at the geometric resonance. At the 
SR frequency (estimated by the asymptotic behavior of the ET dis¬ 
tribution) we underline a weak change in concavity of the curve, in¬ 
dicated by the arrows, (b) PDFs of the escape time, at the frequency 
^ = ^SR = '^/('r)o. (c) PDFs of the ETs at the geometric resonance 
frequency. 


old C: 


Finally, it is possible to design a detection strategy when the 
phase ^ is unknown, that relies on Generalized LRT (GLRT) 
approach, in which a filter bank jointly performs the estima¬ 
tion and detection tuning the initial phase = 2'Kh/H, 
where h = [0, iT — 1]. Thus, it is possible to collect a vector 
of likelihood ratios and the final decision statistics is obtained 
comparing the maximum value of the detectors with a thresh- 


max 

he[0,H-l] 


2N-1 Hi 

Sk{h) < c, 

k=0 Hq 


where the observations are 


Skih) = 0 




(24) 


(25) 
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(again Tg = 2 ttIui). It has been shown that the result¬ 
ing unknown-initial-phase detector for > 10 has nearly 
the same performance as which applies to the coher¬ 
ent (known initial phase) case ||4|, and it is comparable to 
the performances of the noncoherent correlator (standard op¬ 
timum benchmark detector for signals with unknown initial 
phase). On the other hand, the SSC-SRD is computationally 
extremely cheap, in that it only requires binary and/or integer 
arithmetics. 


performances of the detection of Fig. W- at the noise level 
D ~ 0.2 the strobed KC index exhibits a peak in Fig. 

and in the same position the amplitude x, Eq. p7| ), also 
exhibits a peak. A peculiar behavior for the same frequency is 
also observed for the distribution width p, that quantifies the 
degree of synchronization with the input periodic drive ISl: 


IV. DETECTION PEREORMANCES 

From this Section we begin to collect a systematic analy¬ 
sis of the performances of the detection realized through the 
analysis of the ET and the strobed sequence, see Sect. I^The 
numerical simulations have been performed with Euler algo¬ 
rithm taking into account the correction proposed in Ref. Il49l . 
The convergence of the results has been checked both de¬ 
creasing the integration step size and increasing the simulation 
time. Typical results have been obtained with very long sim¬ 
ulations (number of trials for PDF estimation is ^ 10®). The 
observation time (after a transient that has been discarded) for 
signal detection has been set to (T) ~ 4 • 10®. The running 
state analyzed is independent of the initial conditions (initial 
phase (^0 position), and therefore there is no need to repli¬ 
cate the system. 


(here (r^) and (r) are the ETs moments with the applied si¬ 
nusoidal drive). The index ( |28] l measures deviations from the 
exponential case {p = 1) and it is often taken as a signature 
of SR. However, the distortion of the escape time distributions 
of Fig. 1^ are not fateful for signal detection. The analysis of 
the full distribution of the ETs embedded in Eq. 0. not just 
of some parameters as in Eqs. (|27] l and ( |28| ), does not show 
a nonlinear behavior, see Fig. 1^. Thus a physical stochastic 
resonance, such as a peak in the displacement. Fig. or a 
change in the distribution of the escape times. Fig. 1^, does 
not necessarily imply an improvement of the ETs detection 
properties given by Eq. 0. 

We conclude that in the prototypical overdamped bistable 
system, if one wants to exploit SR for signal detection, the 
details of the detection strategy are essential. The analysis of 
the ETs cannot be improved by the increase of noise, while the 
analysis of the strobed dynamics does show an improvement 
at the SR frequency. 


A. First order (overdamped) system 


For the overdamped system governed by Eq.0 the KC in¬ 
dexes df^xc given by Eqs. (14i and ( [2^ , show 

a remarkably different behavior for the two data acquisition 
strategies, see Eig. It is important to notice that the underly¬ 
ing trajectories to be analyzed are the same; the two detection 
strategies are just two different ways to reduce the data, either 
to a sequence of passage times (ET) or to sample the position 
at constant time intervals (SD). 

From Fig. |^it is also clear that also the best performances 
of the detection strategy based on the ETs are below the per¬ 
formances of the strobing in the considered noise interval (we 
measure the quality of the strategies through the dxc peak 
height, and therefore secondary peaks are less interesting). 
The detection performances of the SD deteriorate when the 
frequency increases, while ET performances remain nearly 
unchanged. As shown in Eig. |^, this is consistent with an 
analysis of the displacement amplitude x implicitly defined 

by El: 


{x{t)) = X sin (wf -I- $) . 


(26) 


The quantity x is approximately given by the Equation: 


2a{x‘^)o 
D^4 + oj^t)1 


(27) 


where {x'^)o is the position variance without sinusoidal drive 
(a = 0). Eigure sheds light on the physical origin the 


B. Second order (underdamped) system 

We continue our analysis of the performances of the de¬ 
tection, and now discuss the case of an underdamped system. 
Two regions are a priori interesting: the SR frequency where 
the relation Eq.0 holds ( w ^ 1 for these parameters), and 
the geometric resonance, or the oscillation frequency at the 
minimum of the potential (w ~ 1 in these normalized units). 
The two frequencies are related to interwell oscillations (here 
corresponding to a; ^ 1) and intrawell regimes (here corre¬ 
sponding to w ~ 1). Typical results are collected in Fig. 
|5}i,b, that is the main result of this paper. Figure indicates 
that there are some striking differences with the analogous be¬ 
havior for the overdamped case of Fig. First, at the lower 
frequency corresponding to the SR resonance, the dxc peak 
of strobing better performs compared to the peak of ET, see 
Fig. |5}i. At higher frequencies corresponding to the geomet¬ 
ric resonance, a peak appears that it is only present in the 
underdamped system, see Fig. This peak indicates that 
the ETs better perform respect to the strobed strategy: at the 
geometric resonance. Fig. [^, the role of the best detection 
strategy is exchanged respect to the SR frequency of Fig. 0. 
Thus depending on the nature of the resonance (stochastic or 
geometric) the detection strategy that exhibits the best per¬ 
formances is changed. We ascribe the behavior to the differ¬ 
ent nature of the two peaks: the SR low frequency peak is a 
synchronization phenomenon between noise and determinis¬ 
tic drive that does not accumulate energy in the oscillations. 
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In contrast, the geometric resonance is a standard resonance 
due to energy storage of the oscillations in the well (as such, 
it cannot be observed in the first order system). Thus, the 
two strategies, SD and ET, are best suited for the stochastic 
and geometric resonance, respectively. The two peaks of the 
strategies are found at different noise level (Fig. |^, though 
all parameters are the same. This is possible inasmuch the in¬ 
dex dKC, Eq.(lO), exploits two data sets: one in presence of 
noise alone CHo) and one in the presence of the signal S{t) 
CHi). A maximum occurs when a strategy enhances the dif¬ 
ference, and hence the maximum depends upon the details of 
the employed strategy. In contrast, if one considers SR mea¬ 
sures such as p or x, Eqs. (26) or (28), only based upon data 
sequences containing the signal (Hi), noise most affects the 
sequences at a specific noise level, likely to be the same for 
any indicator (Fig. 3b). 

The analysis of Fig. of the detection features for differ¬ 
ent frequencies of the external drive are further explored in 
Figs. |6|7| and In the neighborhood of ujsr. Fig- 0 the 
performances behavior as a function of the noise identifies a 
peak for the escape based detection, that we compare with the 
frequency ujsr predicted in Eq.Q. The noise intensity at the 
peak of Fig. collected in the inset of Fig. |^, moves to 
higher noise level when the frequency is increased, as pre¬ 
dicted by Eq.([^ that captures the order of magnitude of the 
position of the best performances. We thus confirm that SR, 
the synchronization of the external drive of noise-induced es¬ 
capes, seems to be the physical origin of the peak; however 
SR theory cannot be used for a detailed prediction of the best 
noise input. We remind that in the SR region the peak of the 
ETs in Fig. is not useful for detection, inasmuch the per¬ 
formances are inferior to the strobed performances. 

In Fig. 1^ we investigate the physics and the detector per¬ 
formances around the geometric resonance. The peak of the 
detection performances moves to higher noise when the drive 
frequency increases for the ETs (Fig. and stays at al¬ 
most the same frequency for the strobed strategy (Fig. 0’) 
-the slight change is due to the weak change of the resonant 
frequency, that for nonlinear systems is affected by noise and 
other parameters ll50l . 


To connect the behavior from low to high frequency (i.e., 
from the SR to the geometric resonances) we plot the perfor¬ 
mances of the ETs based strategy in Figs, [^as a function of 
the noise temperature for different frequencies. It is notice¬ 
able that, at variance with the theoretical behavior of the stan¬ 
dard SR Eq.( 261 and with the performances of the detection of 
overdamped systems (Figs. 01 ,b, and c) the dependence upon 
the external frequency exhibits a maximum around w ~ 0.8. 
Also, in contrast with the monotonic behavior of the relation 
Q, the temperature at which the peak occurs for different fre¬ 
quencies first increases (Fig. [^) and then decreases (Fig. [^). 

We underline that the geometric resonant frequency, at 
which the index dj^c peaks, is sensitive to dissipation 7, for 
the response of the system is nonlinear. However, the actual 
peak frequency stays below unity ED. The frequency depen¬ 
dence is also analyzed through Fig. [^, that displays the phys¬ 
ical average escape time as a function of the deterministic sig¬ 
nal frequency. When the drive frequency is in the neighboring 


of the geometric resonance at w ~ 0.8, the escape times are 
shortest, for the average (r) mostly deviates from the average 
escape time without deterministic signal {t)q. In this condi¬ 
tion the performances of the ET strategy are maximized, see 
Fig. [Sf). 

The PDFs of Figs. display other important physical 
features that affect detection. At the lower frequency, corre¬ 
sponding to SR , due to the phenomenon of resonant activa¬ 
tion, the PDF is almost exponential with a cutoff ifSTl . At the 
geometric resonance , due to the phenomenon of dynamical 
resonant activation, the average escape changes and the dis¬ 
tribution develops large oscillations, see Fig. [^, around the 
exponential behavior of the unperturbed (a = 0) system. The 
oscillations reveal that the interwell transitions occur at a spe¬ 
cific value of the signal phase (as observed, for instance, with 
the incoherent acquisition strategy and initial conditions reset 
in Ref. ||50l). This is analogous to the oscillations found in 
experiments on overdamped dynamics of tunnel diodes ll52l . 

In the case of a weak sinusoidal drive, that is the interesting 
limit for signal processing, the signal only introduces a small 
deviation from the unperturbed case instead of a secondary 
minimum ll22l . Actually, an inflection point, indicated by an 
arrow, appears in Fig. [^. 

The decision test, based on Eqs. (13 1^, exploits the os¬ 
cillations of the PDF to distinguish between the perturbed and 
the unperturbed PDFs, and it is therefore most effective at the 
geometric resonance. 

Finally, in Fig. [^we investigate the role of dissipation. 
Figure [TO^ confirms that the detection thorough the escape 
times is favored at low dissipation - indeed for overdamped 
systems it is passed by the strobed strategy. In contrast, the 
strobed based strategy at the lower frequency is little affected 
by dissipation, see Fig. 10 1 . 


V. CONCLUSIONS 

Stochastic resonance, that is an enhancement of the coher¬ 
ent response of a system in spite of an increase of the noise, 
is viewed through the performances of the detection analysis. 
In this framework it is known that, adopting optimal detec¬ 
tion strategies that exploit the full trajectory (i.e.,the matched 
filter), the detection can only improve when noise decreases 
. If instead one is forced to some suboptimal analysis on re¬ 
duced data, as escape times or trajectory discretization (strob¬ 
ing), it has been shown that a trustworthy improvement of the 
performances could occur even in spite of an increase of the 
disturbing noise lUMll- 

This is the premise to investigate the interaction between 
stochastic resonance in the signal detection acception and the 
physical properties of the system. We first construct the like¬ 
lihood ratio test schemes that can be applied to the escape 
times and the strobed discretization. (We note that escapes 
and strobing are, to the best of our knowledge, the only dis¬ 
crete data sample strategies considered in the literature.) The 
measured performances show that: 

1) For overdamped and underdamped systems, around the 
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FIG. 10: Underdamped system, (a) escape and (b) strobed signals for 
two values of the dissipation 7. The escapes peak strongly depends 
upon the dissipation 7. Instead, the strobed strategy is little affected 
by dissipation, as expected for a first order type detection. 


stochastic resonance frequency of Eq. strobing 

overperforms the escape times based strategy 


2) For underdamped systems, at the geometric resonance 
frequency, the escape times based strategy overper¬ 
forms strobing; 


The results reported in 1) are physically consistent with the 
synchronization of the periodic drive and the noise assisted 
leaps over the energy barrier. When the drive period and the 
noise induced jumps are comparable, the escapes distribution 
just remain exponential, see Fig. |^. The escape analysis is 
therefore less effective, as shown in Figs. m- Moreover, this 
type of synchronization can occur for both overdamped and 
underdamped systems, as shown in Figs. 0111, confirming that 
stochastic resonance is inherently a first order phenomenon, 
also in inertial systems. 

The behavior at the geometric resonance, point 2), has also 
a physical origin: it corresponds to the possibility to accumu¬ 
late energy in an eigenmode of the potential. When this con¬ 
dition occurs, i.e. when the periodic drive frequency matches 
the natural frequency of the potential, the distribution of the 
passage times exhibits clear oscillations at the drive period, 
see Fig.J^, and the analysis of the escapes is most effective, 
see Fig. 1^. We remark that the abovementioned detection op¬ 
portunity only arises in underdamped systems that can accu¬ 
mulate energy: The performances of the escapes based strat¬ 
egy improve when dissipation decreases, see Fig. 10 1 , while 



D - Noise Amplitude 

FIG. 11: Threshold setting procedure for the parameter 9 of Fig. [T] 
(a) estimated stochastic resonance frequency [Eq. 0] as a function 
of the threshold 9. (b) closeup ofo.5 < 9 < 0.9 (filled symbols) 
compared with the estimate of losr from the asymptotic exponen¬ 
tial behaviuour of the simulated ETs probability distribution (open 
circles). 


the strobed based detection does not, see Fig. 10 1 . Shortly, 
the physical properties (dissipation, characteristic frequency) 
and the detection strategy are not independent, but deeply in¬ 
tertwined. 


One can conjecture that the prescriptions we have found for 
the analysis of the reduced data can be extended to other sys¬ 
tems as washboard potential m or piecewise linear barriers 
m, that might prove convenient for analytic treatment. As 
a last remark, let us recall that the above results are based on 
a prototypical bistable quartic potential. Results on Josephson 
junctions 11211 l50l and Fabry-Perot pendular interferometers 
G1 are consistent with the above interpretation. 
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Appendix - Numerical procedure to retrieve the Escape Time 
sequence 

It is well known that to determine the passage across a 
boundary with simulations of stochastic differential equations 
is problematic ll49l . In this Appendix we sketch the method 
employed to retrieve the ETs avoiding spurious fluctuations 
across the separatrix around the unstable point of the potential 
Q, a: = 0. In fact, subject to fluctuations, the representative 
coordinate x{t) might cross several times the separatrix, fast 
oscillating around the maximum of the potential before to slip 
down towards a minimum. To avoid such deceitful passages 
that alters the high frequency part of the escape times spec¬ 
trum, we use an effective threshold value ±9, see Fig. [T] If 
the particle is initially in the left hand side of the potential, 
{x < 0 ), we only count a passage across the threshold +9, 
beyond the separatrix in the descending part of the potential 
m- If instead the system is initially in the right hand side 
of the potential (x > 0), an escape is defined as the passage 
across —9. 

The iteration of the previous rules generates an FT se¬ 
quence that is the starting point to estimate the escape 

time distribution at statistical equilibrium, i.e. in the running 
state. The choice of the threshold to decide if the passage to 
the other basin has occurred affects the statistical properties of 
the escape times sequence. The data are recorded when the 
system has reached the running state; the procedure amounts 
to define an escape time as the time to reach the threshold 9 
(—9) starting from the position x = —9 (x = 9 in the basin 
a: < 0 (a; > 0). 

For the first order system this is analogous to the first pas¬ 
sage time from —9 to 9 (or 9 to —9 for the reverse passage). 


For second order systems the analogy with first passage times 
is incomplete, as the initial condition on the velocity is not re¬ 
set. In fact, for a rigorous mathematical definition of the FTs 
it is important the restart process, i.e. to specify the initial 
conditions after each passage (see Ref. Il56l for the role of 
initial conditions). 

To determine a suitable value for the parameter 9, we ob¬ 
serve that the PDF for long escape times values is asymptot¬ 
ically exponential (as expected in the absence of an applied 
external drive, see Fig. Ep ,c), while for low escape time it de¬ 
viates from the exponential distribution. In accordance with 
these observations, we chose the appropriated 9 value impos¬ 
ing that (jJsR, numerically obtained from Fq. 0 . agrees with 
the asymptotic exponential decay ET). 

In Fig. 11a typical estimate of stochastic resonance fre¬ 
quency, computed for the second order system (|^, is dis¬ 
played. The Figure shows that the frequency of the SR (at 
a fixed noise level D) depends upon the choice of the thresh¬ 
old 9 to discriminate a passage from a basin to the other. In 
particular FigfTT^ shows that for 0 > 0.5 a stable ujsr esti¬ 
mated can be achieved. In the close up showed in FigpTjr the 
estimate is compared with the asymptotic exponential value 
of PDF extracted by the simulated FTs distribution; we note 
that a reasonable agreement is found for 9 ~ 0.6, that is the 
value employed in the simulations. 

To guarantee statistical independence from initial condi¬ 
tions, the results are collected after a transient time, typically 
of Ttrans — lOOO(r), where (r) is defined by Fq.([^. We have 
checked that such time is longer than the typical response re¬ 
laxation time of the system. 

Finally, we observe that results similar to Fig. 11 can be 
found for the first order bistable system. 
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